Coherent macroscopic quantum tunneling in boson-fermion mixtures 
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We show that the cold atom systems of simultaneously trapped Bose- Einstein condensates (BEC's) 
and quantum degenerate fermionic atoms provide promising laboratories for the study of macroscopic 
quantum tunneling. Our theoretical studies reveal that the spatial extent of a small trapped BEC 
immersed in a Fermi sea can tunnel and coherently oscillate between the values of the separated 
and mixed configurations (the phases of the phase separation transition of BEC-fermion systems). 
We evaluate the period, amplitude and dissipation rate for 23 Na and 40 K-atoms and we discuss the 
experimental prospects for observing this phenomenon. 
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The tunneling of a macroscopic (or collective) vari- 
able of a many-body system through a classically forbid- 
den region, macroscopic quantum tunneling (MQT), is a 
phenomenon of fundamental interest [1[ and a recurring 
theme in a variety of fields ranging from nuclear (fission) 
and condensed matter physics (e.g. quantum magnets [2|, 
SQUIDs) to quantum optics (macroscopic Schrodinger 
cat states Q and beyond-standard limit measurements). 
Nevertheless, stringent tests under well-understood and 
controlled conditions remain an experimental challenge. 
Cold atom gases, arguably the cleanest and best under- 
stood mesoscopic systems which, furthermore, offer un- 
precedented control knobs such as the ability to vary the 
inter-particle interactions Q, now provide an intriguing 
candidate laboratory for the study of MQT. 

The first cold atom MQT proposals fj| suggested ob- 
serving the collapse of a trapped dilute gas Bose-Einstein 
condensate (BEC) of mutually attracting bosons. How- 
ever, the experimental results [1,0] were either too sensi- 
tive to particle number to distinguish MQT from classi- 
cal collapse [f| , or the analysis was complicated by more 
complex dynamics (such as 'clumping') Q- Evidence of 
coherence (of the many-body system taking on a linear 
superposition of states that correspond to the macro- 
scopic variable residing on cither side of the barrier) is 
even more difficult to gather. Such coherence would be 
more readily observable in the MQT between long-lived 
states, in which case one could set up a coherent pop- 
ulation oscillation between the many-body states. Such 
long-lived states naturally occur in (zero-temperature) 
first order phase transitions in which the order parame- 
ter, which provides the macroscopic variable, can tunnel 
through the barrier of its Landau-Ginzburg potential. In 
the infinite system limit, the coupling between the two 
states rigorously vanishes, but finite- size cold atom sys- 
tems of moderate particle numbers provide, once again, 
a promising candidate to observe the MQT coherence be- 
tween states of different phases, as we show below. 

An earlier proposal to observe MQT between states in 
which the components of a BEC-mixture arrange them- 
selves differently in space, involved a very low coupling 
on account of the small spatial overlap between the sin- 
gle component densities in the different states [8j. In 



this paper, we propose that MQT can be realized and 
its coherence, perhaps, observed in trapped gas mixtures 
of a single-component fermion system and a BEC. Such 
mixtures are currently created 0] e.g. in the sympa- 
thetic cooling scheme in which the colder BEC cools the 
fermions. The tunneling and coherent oscillations that 
we target would occur between states of the mixed and 
separated phases in the phase separation transition of 
the fermion-BEC mixture [To| . Such transitions could be 
accessed by varying the scattering length of the boson- 
fermion interaction [ill ]. 

We consider Nb atomic bosons confined in a spheri- 
cally symmetric harmonic trap (of frequency lot) inter- 
acting with a much larger system of atomic fermions. 
For simplicity we assume the fermions to occupy an in- 
finite volume. The Hamiltonian of the bosons is de- 
scribed by the standard Gross-Pitaevskii (GP) form [l|, 
i.e., with inter-particle interactions described by a con- 
tact potential (oc Ass<5(r — r')), which we choose to be 
repulsive (Xbb > 0) We assume that the interaction of 
bosons with fermions is also contact-like, contributing 
A_bf| x I / s| 2 | v E'f| 2 to the Hamiltonian density, where Xbf 
is the fermion-boson coupling constant. Furthermore, 
all fermions occupy in the same spin state so that the 
short-range inter-fermion interactions do not contribute 
by virtue of the Pauli exclusion principle. 

We are interested in the dynamics of the reduced sys- 
tem of bosons described by the functional 
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where So is the action of the bosons alone, Sbec — 
J dr{h J d 3 r^' B ^'|j — H B ec): an d the second term is 
a contribution due to the interaction of bosons with 
fermions; \xp is the chemical potential of the fermions. 
Here and throughout the paper we will be utilizing 
the imaginary time (Matsubara) representation, unless 
stated otherwise. An explicit evaluation of the second 
term is a challenging task. However, here we are inter- 
ested in the dynamics of the slow breathing mode of the 
BEC ^ B > which can be treated in the self-similar density 
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approximation. This dynamics describes the longitudinal 
expansions (and contractions) of the condensate. Finite 
size effects such as the appearance of a non- vanishing ex- 
citation energy (gap) can decouple this mode from other 
excitation modes. Hence, ^ B peaks at small frequen- 
cies (d) and small wavevectors (q), giving a ^° B that is 
a slowly varying function of spacial and temporal coor- 
dinates. In such case the Trlog[...] in Eq. Jl} can be 
evaluated within the Thomas-Fermi approximation. A 
straightforward zero-temperature calculation yields 
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where k F is the Fermi wavevector. Eq. ^ represents an 
additional term in the Gross-Pitaevskii (GP) equation, 
SSgp/S^b = 0) resulting from interaction with fermions. 
In order to analyze the physical meaning of Eq. (U) let 
us expand it in powers of ^ B . The first nontrivial con- 
tribution is a term -2A'|* b | 2 *b, A' = (X 2 BF k F /4Tr 2 ^ F ), 
which corresponds to the attraction between bosons me- 
diated by interaction with fermions. For nonzero, but 
small uj and q there is an additional term (of the or- 
der of X BF ) related to the dissipation of the condensate 
due to the Landau damping, as we discuss below. The 
next order yields ?7|*b| 4 *b, V = (k F \ BF /8ir 2 n F ). Un- 
like the previous term this one is positive, and represents 
reduction in the effective boson-boson attraction due to 
depletion of fermions in the regions of high density of 
the bosons. The next order terms (in Xbf) prove to 
be unimportant as can be verified directly from Eq. Ip[). 
Therefore we will replace the potential energy contribu- 
tion in GP equation given by Eq. ^ by the two terms 
discussed above fl2j |. 

To analyze the dynamics of the slow (breathing) mode 
described by the Hamiltonian 
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we apply the time dependent variational principle. Since 
we are interested in ground state properties of Eq. (5) we 
use a spherically symmetric Gaussian trial wavefunction 
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parameterized by a dimensionless parameter x that char- 
acterizes the BEC's spatial width in units of the zero 
point motion amplitude Rq = (h/2niBi^T) 1 ^ 2 ■ Substitu- 
tion of this wavefunction into Eq. (5) yields the following 
dependence of the ground-state energy E° on x: 
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where a = N B (X' - 
4iV|r//(3 5 / 2 7r 3 i?6^ T ) 
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FIG. 1: Density profiles of the bosons (solid lines 1 and 2) and 
corresponding density profiles of the fermions (dashed lines 
1' and 2') in two metastable states: (1) with fermions having 
zero density at the center of the trap ( "separated phase" ) and 
(2) with nonzero density of fermions ("mixed phase"). The 
dotted line represents schematically an effective potential for 
the breathing mode of the bosons. 



a, i.e., for a < a cr = 32(2/5) 1 / 4 /I5 ~ 1.69, E may de- 
velop two competing minima, depending on the value of 
the /3-parameter. The energy barrier separating the min- 
ima is caused by the same effect as the barrier appearing 
in the description of a BEC with attractive interactions: 
it arises due to the competition between the kinetic and 
the interaction energies, i.e, the first and the third terms 
in the right-hand side (rhs) of Eq. (J5J. In the absence of 
the last term in the rhs of Eq. ((5j) the state in this well 
would have been metastable - the energy would tend to 
— oo at x — > 0. The 1/x 6 term stabilizes the system: for 
small x this term rapidly increases, giving rise to another 
minimum of Eq(x), now due to the competition between 
the last two terms in the rhs of Eq. ([5]) . 

At certain values of a and (3 the two minima of E$(x) 
will have the same energy, and the ground state of the 
system becomes degenerate. Since our system is finite, 
this degeneracy will be lifted by the quantum tunnel- 
ing transition between the two states. Such mechanism 
has been suggested to be the dominant decay process for 
condensates with attractive interactions between parti- 
cles Q . The tunneling corresponds to the low energy ex- 
citations of the breathing mode, described by the wave- 
function in Eq. It has been shown in [5( that by 
accounting for the superfluid motion of the condensate 
(which can be done by introducing a phase-factor e 1 ^ for 
the wavefunction in Eq. (j4|) and requiring the superfluid 
velocity v s = (h/m B )V<fi to satisfy the continuity equa- 
tion) one obtains an effective action for the breathing 
mode of the condensate 
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where E (x) is given by Eq. §5§ and m = 3m B N B R,Q/2. 
Thus the dynamics of the ground state wavefunction of 
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the condensate is that of a quantum particle of mass mo 
moving in the potential E (x). 

A direct analysis of the Shrodinger equation corre- 
sponding to Eq. ([S]), however, is quite cumbersome since 
the two wells are generally quite asymmetric. Instead 
we choose an alterative route: we compute the ground 
state energy and obtain the tunneling rate by numerically 
solving the time-independent GP equation, SH/S^b = 
E^b, where H is given by Eq. ([3]). The latter approach 
also serves as an independent justification of the varia- 
tional method and confirms that macroscopic quantum 
tunneling, QMT, is the mechanism that causes the tran- 
sition between the two states of the condensate. Upon 
substitution ^>b ~ 4>/r, the time-independent GP equa- 
tion can be cast in the form 
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where the </>(a;)-function is normalized to unity, a = 
( 7 r/2) 1 / 2 a, b = 3 5 / 2 tt/3/16, and x = r/R , fi = E/fuv T . 
We find the ground state numerically by replacing the 
rhs of Eq. (J7J) by ~d T cj) and propagating </> in the imag- 
inary time t until it converges to the ground state <^>o 
(or *& B ). We then evaluate the ground state energy ac- 
cording to Eq. ([3]) and present the results in Fig. 2(a) 
as a function of the 6-parameter for different values of a. 
Fig. 2(b) shows the dispersion of the ground state width, 
(1/iVfl) / d 3 r|\f r 'g(r)| 2 r 2 , as a function of those same pa- 
rameters. For a < a cr — 1.83 the ground state energy 
and dispersion undergo a sharp crossover between the 
state with compressed and expanded BEC wavefunctions 
(corresponding to the phase separated and mixed states) 
as functions of b. Note that the value a cr corresponds 
to the value of a cr = 1.46, which is quite close to the 
above critical value of 1.69 obtained from the variational 
approach. The dependence of ground state energy near 
the critical value of a is shown in the inset of Fig. 2(a). 
Clearly the ground state energy exhibits avoided level 
crossing, which is in accordance with the above conjec- 
ture (e.g., Eq. ([6])) of macroscopic quantum tunneling 
between the two local energy minima. 

The value of the tunneling matrix element A between 
two local "ground" states t\ and £2 can be deduced by 
fitting the calculated energy curves in Fig. 1 with the 
standard expression [l3j |. e = (ei + ea)/2 — [(ei — £2) 2 /4 — 
A 2 ] 1 / 2 and assuming that in the vicinity of the point 
of crossover both e\ and e 2 are linear functions of pa- 
rameter b. For a = 1.81 one finds A ~ 1CP 4 x Hujt, 
while for a = 1.82, A ~ 10~ 2 x Hu>t- Assuming that 
the ground state wavefunctions have Gaussian shape, 
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from Fig. 2 one finds 

R 2 )/2 ~ 0.85i? for both a = 1.81 and 
a = 1.82, and 5R = \Ri - R 2 \ ~ 0.09i? for a = 1.81 
and 5R ~ 0.03i?o f° r a = 1-82. For a typical value 
of the trapping frequency vt = lQ 2 Hz (wt = 27t^t), 
the two tunneling rates are Ai.gi/fi, = 10~ 2 s~ l and 
Ai.82/^ = 10 2 s _1 . Since the value of R for most 
trapped atomic BEC's is of the order of a few microns, 
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FIG. 2: (a) Dependence of the ground state energy of the BEC 
(per particle, in units of Hlot) as a function of parameters a 
and 6; (b) Dispersion of the ground state spatial extent as a 
function of the same parameters. 



the difference between the radii of the two condensate 
states 5R is submicron. Such small variation may be 
difficult to observe in situ by optical means. However, 
the expansion process that takes place in time-of-flight 
measurements after the trap potential is shut off and the 
expanding atoms are observed, has successfully magnified 
small distance features in other experiments. 

Role of dissipation: The above analysis determines 
the tunneling rate, but does not address the question 
whether the tunneling process is quantum coherent. Will 
the probability of the system to occupy one of the two 
macroscopic states oscillate in time as cos 2 (Ai/fi)? The 
fermions not only provide the BEC with the effective in- 
teraction, they also cause fluctuations which can destroy 
the macroscopic quantum coherence. To evaluate the 
effect of fluctuations, it is sufficient to consider the first 
non- vanishing frequency-dependent contribution into the 
effective action of the bosons coming from the perturba- 
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tive expansion of the Trlog [...] term in Eq. (fTj): 

Here ps(q, lo) is the Fourier transform of ps(r, t) and 
Xo is the response function of the non-interacting 
fcrmions. In the small frequency domain xa = 
(l/4ir)[H 2 k F /(nfi F ) + m 2 F \ui\/(h 2 q)]. The frequency- 
independent part of Xo has already been incorporated 
in the effective interaction between bosons, i.e., A'^sj 4 
term in Eq. ([3]). The second term in \o is responsible 
for damping. To quantify its role we employ a two- 
state approximation in describing the tunneling dynam- 
ics. In this representation the tunneling is described by 
the Hamiltonian H tun — Aa x , where a x is a Pauli ma- 
trix with non-zero off-diagonal elements, and the posi- 
tion operator, i.e. the spatial width of the ground-state 
BEC wavefunction, is given by R = R + (SR/2)a z , where 
a z is the diagonal Pauli matrix (with ±1 along the di- 
agonal). The dissipative part of the action for H tun 
can be derived from Eq. © by substituting a Gaus- 
sian ansatz, p B _(r,t) = N B /{Tr 3 / 2 R 3 (t)) exp {-r 2 /R 2 (£)], 
where R(t)_= R+ (5R/2)<r z (t), a z = ±1, into Eq. ©. 
For SR <C R one obtains 

S diss = jh J dTdT'a z (r)<j z (r')(T - r')-\ (9) 

where 7 = N 2 \ 2 BF m 2 F 5R 2 /[2(2irh) A R i }. Eq. ©, to- 
gether with H tun defined above, describes dissipative dy- 
namics of a two-state system. Such dynamics has been 
extensively studied in connection with macroscopic quan- 
tum tunneling of a superconducting phase in Josephson 
junctions, and is known to depend critically on the value 



the parameter 7. Specifically, for 7 > 1 the two-state os- 
cillation is always overdamped and at zero temperature 
it exhibits localization as a result of quantum fluctua- 
tions [HI • It is therefore instructive to evaluate 7 for our 
situation. For estimates we consider an atomic mixture 
of 23 Na (bosons) and 40 K (fermions), which have natural 
scattering lengthes clbb — lnm (Xbb = 4,Trh 2 aBB /tub) 
and asF — 4nm (Xbf — 2-Kh 2 aBF[{]-/mB) + (I/wf)])- 
For these data we obtain a critical value of N B r ~ 12400 
(again for — l0 2 Hz) and the fermion density n° F ~ 
7.4 x 10 15 cm~ 3 . Then, for a — 1.81 we obtain 7 181 ~ 1.1, 
which corresponds to the localized case (at T = 0), 
whereas for a — 1.82 one gets 71.32 — 0.1. In the high 
temperature limit (for ksT > A) the relaxation rate T 
can be expressed in terms of 7 as fir = ix^ksT [3], 
and therefore coherent (underdamped) oscillations can 
be observed for T < Ai. 82/(71.82^5) = 0.5nK. The 
situation can be improved, however, if one utilizes a Fes- 
hbach resonance [4| to increase the clef scattering length. 
For example, for clef = 80rtm one finds N B r ~ 25 and 
~ 2.6 x 10 n cm~ 3 , and 7l . 82 ~ 2.5 x 10" 4 . For 
such parameters coherent oscillations can be observed for 
T -^i 0.2fiK, which is easily observable. A low particle 
number also reduces the uncertainty of an atomic count- 
ing measurement that can be carried out in the time-of- 
flight procedure [HI]. 

In summary we argue that a trapped boson-fermion 
mixture can exhibit MQT tunneling and coherent oscil- 
lations. Our studies indicate that MQT can be observed 
in 23 Na and 40 K atomic mixtures of sufficiently low tem- 
peratures. 
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